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"  —Thie  paper  deals  with  the  development  of  a  method  for 
calculating  supersonic  flow  past  blunt  bodies  with  a  detached  shock 
wave,  which  was  proposed  first  by  G.  F.  Telenin  of  the  Scientific  and 
Research  Institute  of  Mechanics  of  Moscow  State  University  in  1961. 

The  first  section  contains  a  formulation  of  the  problem  of  supersonic 
flow  past  a  blunt  body  of  an  arbitrary  shape  and  an  Investigation  of 
basic  gas-dynamics  equations.  A  method  is  given  for  the  selection  of 
a  differential  system  for  solving  a  boundary  problem.  This  is 
followed  by  sections  on  the  construction  of  an  algorithm  for  the 
numerical  solution  of  a  boundary  value  problem  in  a  mixed  region  and 
on  the  method  of  solving  a  nonlinear  problem.  The  final  section 
gives  examples  Illustrating  the  solution  of  different  problems  of  flow 
past  bodies  of  different  shapes  with  a  detached  shock  wave.(  Orig.  art. 
has:  35  citations  in  the  references  at  the  end  of  the  article, 

16  figures  and  72  equations. 
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APPLIED  PROBLEMS 


METHOD  FOR  CALCULATION  OF  SUPERSONIC  FLOW  PAST 
BLUNT  BODIES  WITH  A  DETACHED  SHOCK  WAVE 


S.  M.  Silinskiy,  G.  F.  Telenin 
and  G.  P.  Tinyakov 

In  1961  at  Moscow  State  University  G.  F.  Telenin  proposed 
a  method  for  calculating  supersonic  flow  past  blunt  bodies  with  a 
detached  shock  wave  [23,  24].  During  1962  and  1963  systematic 
investigations  were  conducted  on  an  ideal  gas  flowing  past  ellipsoids 
with  a  different  ratio  of  semiaxes  and  bodies  with  an  analytic 
contour,  close  in  form  to  a  face  with  a  small  degree  of  rounding 
of  angles  [27,  32],  bodies  with  a  concave  contour  in  the  region 
of  critical  point,  bodies  with  discontinuity  of  contour  curvature 
in  a  subsonic  unit  and  contour  fracture  in  the  sonic  point  (type 
of  face),  and  flat  bodies  [32,  33];  flow  past  bodies  with  a 
various  form  with  M  number  striving  for  infinity  and  the  aaiabatic 
index  striving  for  unity  [27];  flow  past  a  sphere  by  equilibrium 
carbon  dioxide  [26]  and  past  bodies  of  different  form  by  equilibrium 
air  [25,  34];  flow  past  a  sphere  by  nonequilibrium  air  [28], 
taking  into  account  kinetics  of  excitation  of  fluctuations  in  oxygen 
and  nitrogen,  the  flow  of  dissociation  and  exchange  reactions 
(including  0^  +  **  2N0)  and  ionization  with  the  help  of  atomic 

collisions  (N  +  0  ■  N0+  +  e,  N  +  N  ■  N^  +  e,  0  +  0  ■  0^  +  e),  satis¬ 
factorily  describing  the  process  up  to  M  'v  30  [29  ,  30],  and  spatial 
flowing  around  with  a  detached  shock  wave  [2,  35].  At  the  same  time 
work  cr.  investigation  of  the  method  was  conducted. 

In  the  first,  second,  and  third  sections  of  this  article  the 
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problem  of  supersonic  flow  past  bodies  with  a  detached  shock  wave 
is  formulated,  and  on  the  basis  of  an  analysis  of  properties  of 
solutions  of  model  linear  equations  of  elliptic  and  mixed  types 
the  basic  ideas  of  the  method  are  expounded.  In  the  fourth  section 
a  numerical  method  for  solving  nonlinear  problems  is  described. 

In  the  fifth  section  examples  are  given  which  illustrate  the  applica¬ 
tion  of  the  method  for  solving  different  problems  of  flowing  around 
with  a  detached  shock  wave 

§  1 .  Formulation  of  Problem 

We  will  examine  the  flow  past  an  axisymmetrical  body  by  a 
stream  of  an  ideal  gas  with  the  parameters  Vfc,  pm,  p^. 

Using  the  stream  function  \p ,  the  system  of  equations  for 
gas  dynamics  in  a  spherical  system  of  coordinates  r,  6,  <J)  can 
be  presented  in  the  following  form: 


Here  '  '  _)  .  1  >j  1  j,  a,  projection  of  velocity  vector  W  on  unit 
vectors  j  and  relative  to  V  ;  p,  p  -  pressure  and  density 

A  p  V  I  I  Id  A 

relative  to  p  v  and  p  respectively;  r  is  relative  to  the  character- 
®  max  00  j  2 

istic  dimension  L,  and  ^  to  p  V  L  . 

*  ®  max 

The  basic  difficulty  is  in  the  transonic  nature  of  the 
problem,  since  the  solution  has  to  be  looked  for  simultaneously 
In  the  subsonic  and  in  a  certain  part  of  supersonic  regions.  As 
is  known,  Lr.  the  transonic  region  the  V  ~ :  c  -.:r.  ‘-quuti  ..n  ;  s.  souna 


-•  j 


and  is  the  most  studied  model  of  gas  dynamics  equations  in  the 
hodograph  plane.  It  has  been  used  to  help  in  studying  the  formulation 
of  ma  .y  transonic  problems.  It  is  natural  to  lean  on  these  results 
also  in  examining  the  problem  of  flow  past  with  a  detached  shock 
wave.  According  to  the  results  obtained  in  the  present  work  (Pig. 
la),  a  geometric  picture  is  given  for  flow  past  a  sphere  by  an 
ideal  gas  (y  «=  1.4)  at  M  ■  6.  The  figure  shows  the  shock  wave 
ABC,  sonic  line  DBQ,  maximum  characteristics  DCK,  limiting  fclje 
region  of  influence  of  the  supersonio  fraction  of  flow  behind 
the  shock  wave  on  subsonic.  Also  shown  are  a  series  of  characteristics 

i 

BK,  MN,  etc.,  of  the  first  family,  emanating  from  the  sonic  line* 
and  the  limiting  line  QP,  starting  at  the  point  of  intersection  of^ 
sonic  line  with  the  axis  and  representing  the  envelope  of  character¬ 
istics  of  the  first  family.  Solution  for  the  region  located  before 
the  shock  wave  is  obtained  numerically  by  continuation  from  data 
behind  the  shock  wave  with  the  help  of  the  same  program.  With 
a  decrease  of  M  number  and  an  increase  of  body  flatness  the  form 
and  arrangement  of  mutual  disposition  of  the  compression  shock, 
sonic  line,  and  maximum  characteristics  are  changed.  For  clarification 
of  the  main  positions  we  will  limit  ourselves  to  an  examination 
of  the  simplest  case  which  is  depicted  in  Fig.  la.  In  Figure  lb 

a  representation  of  flow  in  the  hodograph  plane  v  ,  v  (corresponding 

x  y 

points  are  designated  by  the  same  letters  with  primes)  Is  given 
for  this  case. 


Fig.  la.  Fig.  lb. 
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Solution  of  system  (1.1)  should  be  subordinated  to  correlations 
for  compression  shocteT 


>J  •'III  o  —  0t  COS  0,  0  -*  Uy  CO 3  0  -J-  ux  o!'t  l), 

of !-  vV)i 

V)  v,  <  I.  Y  +  1  V  cr  y  J 


Or  & 
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^  jmm _ ^ 


y  —  I  g*  vt 

m«  .  rm 

I - ; —  cos*  u 

r* 

im* 

o  ~  -  - .  —-•**—  A.  siu1 0. 

Y  o  V  c 

*  ri.  .* 


(1.2) 


Here  v  and  v  are  projections  of  velocity  vector  W  on  rectangular 
x  y 

axes  x  and  y  (Pig.  la),  and  r  ■  r  (0)  and  a  the  equation  of  shock 
wave  and  angle  of  inclination  of  it  to  the  axis  of  symmetry  x, 
to  conditions  of  symmetry  on  segment  of  axis  AE  and  the  boundary 
conditon  on  body  contour 


rT 

U  s»  0— — i 
rt 


(1.3) 


where  r  «  rT(0)  -  equation  for  body  generator. 

Prom  the  relationships  (1.2)  it  is  clear  that  if  the  form 
of  shock  wave  is  given,  then  u(i')  is  a  known  function. 


In  regions  A'B’K'D’E'A’  in  the  hodograph  plane  a  boundary 
value  problem  of  the  Frankl-Tricomi  [3]  type  appears:  boundary 
conditions  are  assigned  on  the  elliptic  segment  of  the  boundary 
of  region  B'A’E’D*  and  on  the  noncharacteristic  curve  B 1 C ' ,  located 
inside  the  characteristic  triangle  B’K'D’  which  is  adjacent  to 
the  parabolic  line  D'B*.  Line  A'B'C',  which  is  a  form  of  compression 
shock,  Intersects  each  of  the  characteristics  only  once  inside 
triangle  B'K'D'.  Line  E'D',  a  form  of  body  contour,  is  unknown 
beforehand  and  should  be  determined  during  the  solving  process. 

In  an  analogy  with  the  Tricoml  equation  and  flat  potential  rrcrlems 
of  gas  dynamics  one  may  assume  that  in  the  Pedograph  plane  the 


correlations  for  section  A'C'  o'  compression  shock,  conditions 
of  symmetry  cn  the  segment  of  axis  A  *E  • ,  and  the  boundary  conditions 
on  part  of  body  contour  E’D’  uMquely  determine  the  solution  in 
region  A’B’K’D’E'A’ .  However,  a  physically  possible  flow  does 
not  satisfy  every  solution  in  &  hodograph  plane,  since  mapping 
from  a  hodograph  plane  into  a  physical  plane  can  be  single-valued, 
which  leads  to  peculiarities  in  tne  physical  plane.  In  the  subsonic 
part  of  the  hodograph  plane  mapping  into  a  physical  plane  does 
not  have  any  peculiarities  beside4?  the  isolated  branch  points. 

In  the  supersonic  region  limiting  lines  can  appear,  along  which 
the  fastening  of  two  sheets  of  the  physical  plane  occurs.  Along 
limiting  lines  the  Jacobian 


D  , 

^  K°v> 

beconjes  zero,  or  which  in  our  easy  is  equivalent,  derivatives  of 
velocity  components  become  infinite  . 

In  order  that  the  solution  in  the  hodograph  plane  for  region 
A'B'K’D'E'A'  in  the  physical  plane  ii:i  the  >egiori  AftKDEA  is  satisfied 
by  the  real  flow,  it  is  necessary  to  place  an  additional  condition 
for  the  desired  solution:  nonrevevsni  of  the  Jacobian  D  to  zero, 
or  what  is  the  same,  boundedness  of  derivatives  of  vectorial  velocity 
components  in  the  triangular  region  HKD.  In  accordance  with  this 
the  condition  of  uniqueness  in  reference  to  the  problem  of  supersonic 
flow  past  a  blunt  body  with  a  detached  shock  wave  can  be  formulated 
in  the  following  vay:  in  the  clast  of  functions  with  a  continuous 
limited  derivative  of  relationship  foe  section  AC  of  the  compression 
shock,  the  conditions  of  symmetry  for  segment  of  axis  AE.  and  the 
boundary  condition  for  a  section  of  .,ody  cento*;?*  uniquely  determine 
the  solution  of  a  system  of  differential  equations  of  gas  dynamics 
in  the  region  ABODE  of  the  physical  p.l4ne  (see  Fig.  1;. 

It  is  necessary  to  rote  that  sir.ee  fulfillment  of  the  character¬ 
istic  equation  along  the  characteristic  curve  is  a  corollary  of 
differential  equations,  then  assignment  of  this  equation  along 
limiting  characteristic  DK  or  some  other  characteristic  can  produce 
nothing  for  separation  of  the  class  of  solutions  with  a  limited 
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derivative  (such  a  condition  is  D  /  0  in  the  supersonic  region), 
and  consequently  also  for  formulation  of  conditions  of  uniqueness. 

The  appearance  of  a  limiting  line  is  conditioned  by  the  boundary 
conditions  of  the  problem.  The  characteristic  equation  contains 
derivatives  only  along  the  curve  and  is  fulfilled  for  any  characteristic 
independent  of  whether  or  not  it  has  a  common  point  with  the  limiting 
line  and,  consequently,  whether  or  not  the  derivatives  which  are 
normal  to  characteristic  remain  limited  (for  example,  formation 
of  a  limiting  line  during  supersonic  flow  around  a  concave  angle). 


Making  further  use  of  well-known  theorems  for  hyperbolic  equations, 
our  result  is  that  in  the  class  of  functions  with  a  limited  derivative 
the  solution  is  determined  uniquely  in  a  wider  region  ACGRFDE, 
if  the  conditions  for  compression  shock  and  for  the  body  are  assigned 
accordingly  up  to  points  G  and  F  (Fig.  2). 


If  the  body  contour  up  to  point  T  is  an  analytic  curve,  t her., 
as  it  follows  in  the  region  which  is  limited  from  above  by  the 
characteristic  emanating  from  point  T  the  field  of  gas-dynamic 
parameters  behind  the  shock  wave  not  only  satisfies  conditions 
of  uniqueness  but  also  is  analytic.  The  angular  point  and 
discontinuity  of  curvature  (or  leading  derivatives)  on  the  contour 
in  the  subsonic  sector  of  flow  lead  to  appearance  of  isolated  singular 
points  on  the  surface  of  the  body. 


Let  us  trace  in  region  ACGRFDE  an  analytic  curve  OS?  which 
does  touch  upon  a  characteristics  in  one  point.  Then,  if  the  unknown 
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solution  is  analytical  In  a  sufficiently  large  region,  then  in 
the  region  FSGTF  which  is  bounded  above  by  the  arbitrary  curve 
GTF,  there  exists  an  unique  solution  of  the  Cauchy  problem  with 
initial  data  for  line  GSF  and,  consequently,  in  the  class  of  analytic 
functions  the  above-mentioned  conditions  of  the  boundary  value 
problem  determine  the  unique  solution  in  region  ACGTFDE. 

§  2 .  Selection  of  a  Differential  System  For  Solution 
of  the  Formulated  Boundary  Problem 

For  this  purpose  it  is  possible  to  use  the  method  of  establish¬ 
ment  or  some  iteration  method  for  solving  a  system  of  difference 
equations  approximating  a  system  of  stationary  equations  (1.) 
simultaneously  in  all  the  subsonic  and  transonic  regions  (i.e., 
the  solving  of  nonstationary  equations).  However  both  these  methods 
require  a  great  volume  of  computer  storage  and  the  carrying  out 
of  a  large  volume  of  calculations,  since  they  essentially  reduce 
the  problem  to  a  three-dimensional  one.  This  limits  the  selection 
of  high  speed  machines  for  the  resolution  of  problems  of  supersonic 
flow  past  bodies  with  detached  shock  wave.  Therefore  it  is  natural 
to  attempt  [4,  6,  7]  to  construct  a  numerical  algorithm,  using 
the  resolution  of  the  Cauchy  problem  in  some  form  or  other.  The 
difficulty  is  that  in  a  general  case  the  Cauchy  problem  is  tactless 
for  an  elliptic  region,  i.e.,  solutions  exist  with  Cauchy  data 
for  a  length  of  boundary  of  the  region  differing  as  little  as  desired 
but  which  differ  as  strongly  as  desired  Inside  the  region. 

However,  in  a  number  of  works  (see,  for  example,  [8,  9]  it 
is  shown  that  the  problem  becomes  correct  in  a  limited  region  if 
one  were  to  narrow  down  the  class  of  solutions  examined.  In  the 
class  of  bounded  functions  the  problem  is  correct  In  the  sense 
of  the  root-mean-square  norm,  and  if  the  function  is  moreover  analyti 
then  in  the  usual  sense. 

Considering  the  analytic  nature  of  the  desired  solution,  it  is 
natural  to  expect  that  by  selecting  an  approximating  system  which 
takes  this  analyticity  into  consideration  to  a  sufficient  degree, 
it  is  possible  to  obtain  a  convergent  and  sufficiently  stable  method 
for  the  numerical  resolution  of  the  problem. 


We  will  consider  the  selection  of  an  approximating  system  in 
an  example  of  the  Laplace  equation: 


di*' 


►  . 
I 


(2.1) 


Let  it  be  necessary  to  find  the  solution  of  a  Cauchy  problem  in  the 
region  ^  ■  ’•  with  initial  data  (Fig.  3) 

•f(*.  0)  •»«(«).  v(a).  (2.2) 

Let  us  assume  that  in  a  certain  wider  region  ABCD  the  desired 
solution  <Mx,  y)  is  a  harmonic  function  with  a  limited  gradient 
|  grad  <p  |  ±  M,and  the  initial  data  (2.2)  on  segment  KL  is  less  than 
( 1/2 )m. 


sums  m 
n^cm 


Fig.  3. 

Then  for  the  function  of  the  complex  variable  f(z)  *  t  ( x ,  y)  + 
+  lt^(x,  y)  ( ( 0 . 0 )  *  0)  we  have  the  estimate  .'i  i.  a'  ./  > 

in  the  region  ABCD  (e  <  H)  and  !  u)  U  for  the  initial 

line  KL.  Using  the  Corleman  lemma  [8,  9]  we  obtain  tne  estimate 
for  resolution  of  the  Cauchy  problem  in  the  region  I  a  J 
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Formula  (2.3)  -hows  that  the  gr-ater  a  and  ii  -  b  are,  (tin; 
further  peculiarities  are  located)  then  the  more  rapidly  decreases 
the  solution  of  4>(x,  y)  in  the  entire  examined  ( crosshatched )  region 
with  a  decrease  of  m.  The  value  E  /  0  is  not  essential. 

Now  let  us  consider  on  the  same  nodel  the  question  of  the 
selection  of  an  algorithm  for  the  numerical  resolution  of  the  Cauchy 
problem  in  an  elliptic  region. 

If  one  were  to  conti  ue  the  entire  solution  $(x,  y)  = 

=  4>  ( x ,  n,  y)  +  i<j>2(x,  n,  y)  in  the  complex  region  C  *  x  +  iy, 
then  at  a  fixed  x  (x  -  parameter)  determination  of  4^  and  is 
reduced  to  resolution  of  the  Cauchy  problem  in  the  plane  n,  y  for 
a  hyperbolic  equation  (in  the  examined  case  wave  [10,  11])  with 
initial  data 

•fi ( * .  >i)  »Ti ('.  n).  «fi, (*.  >i)  ~  M*.  n).  ^(x,  >])  ir,(.o  n), 

T. >,  (  <■.  »l)  *t,(r, \\). 


Here  •  ij)  ;  »'*  •))  and  x(Q  -  t,  (x,  rtf  »i).  are  analytic 

continuations  of  Initial  data  into  the  complex  regions  c  *  x  +  in. 

Thus,  In  the  three-dimensional  region  x,  n,  y  it  is  possible 
to  construct  a  stable  algorithm  for  the  numerical  resolution  of 
the  Cauchy  problem  and  to  find  the  solution  of  the  initial  problem 
<f>(x,  n)  *  ^(x,  0,  y) . 

Such  an  approach  has  two  deficiencies  from  the  point  or  view 
of  practical  use  in  the  numerical  resolution  of  a  Cauchy  problem 
in  an  elliptic  region.  First,  the  process  of  continuation  of 
initial  data  into  a  complex  region  is  reduced  to  resolution  of 
the  Cauchy  problem  for  the  Laplace  equation.  For  example,  for 
determination  of  o^(x,  n)  it  is  necessary  to  solve  the  problem 
with  Initial  data 


•Y.('.O) 


d,  (  V.  0)  O  (  <), 


Just  as  in  a  general  case,  the  problem  is  reduced  to  only  a  simple 
Cauchy  problem  in  an  elliptic  region,  and  in  the  examined  case 
of  a  model  equation  this  problem  is  equivalent  to  the  initial  one. 
Secondly,  there  is  a  strong  increase  in  the  volume  of  calculations 
since  actually  again  a  problem  is  being  solved  in  the  three- 
dimensional  region  x,  n,  y.  However,  both  these  difficulties  can 
be  surmounted  by  the  selection  of  a  suitable  difference  system 
in  the  region  of  "physical"  variables  x,  y,  giving  sufficient 
consideration  to  the  analytic  nature  of  the  desired  solution. 


The  general 
Laplace  equation 


solution  of  the  Cauchy  problem  (2.2)  for  a 
(2.1)  can  be  presented  in  the  form 


i ) 


m,v)  -  (>04 

a 


where  o(z),  T(z)  -  analytic  continuations  in  the  complex  region 
z  =  x  +  iy  of  functions  o(x)  and  T(x). 


Let  us  assume  that  for  the  initial  segment  tr.e  r  odes 

2n  +  1  are  selected.  Based  on  the  values  of  functions  in  one  r.r-'t 
for  tr.e  initial  data  we  construct  Lagrangian  interpolation 
polynomials  in  the  complex  region: 


I  ’  it 


where  a,  and  b  are  linear  functions  of  o(x.)  and  t  ( x  .  /  value:- 

J  » J  j 

functions  in  the  nodes.  If  functions  o(x)  and  t(x;  are  int , 
then  during  arbitrary  location  of  nodes  the  sequence  of  Inter¬ 
polation  polynomials  (2.5)  converges  evenly  in  any  finite  regie, 
of  the  complex  variable. 


If  functions  o(r, )  and  t(c) 
plane  and  Interpolation  points 
polynomial,  tro.-  tr.<  interpolat 


have  peculiarit 
coincide  with  ze 
in.a  proce.;. 


1  o  :  In  the 
roe.'  '/'  t he¬ 


ron  ve 
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in  any  closed  region  lying  inside  the  ellipse  r  - 1 1  p*. 

Tf  all  nodes  coincide  with  the  principle  of  coordinates,  then 
polynomials  (2.5)  convert  into  segments  of  power  series,  and  the 
region  of  uniform  convergence  is  any  closed  region  which  is 
internal  with  respect  to  circle  R  =  | £ |  =  R* .  Ellipse  p  =  p*  and 
circle  R  =  R*  pass  through  the  nearest  singular  point  (Fig.  *0 . 1 


Fig.  4 


GRAPHIC  NOT 
REPilCulCIBLE 


Let  us  note  that  the  use  of  nodes  which  are  disposed  on  the 
segment  expands  the  convergence  region  of  the  interpolation 
process  in  direction  of  the  initial  line  without  changing  its 
dimension  in  the  direction  n,  when  the  peculiarity  is  disposed 
at  x  =  0  and  a  fixed  n.  Inside  the  region  of  uniform  convergence 
the  interpolation  process  gives  an  effective  stable  (^n  the  examined 
class  of  functions)  algorithm  of  continuation  o'  initial  data  in 
the  complex  region. 


Placing  (2.5)  in  (2.4),  we  obtain  an  approximate  solution 
of  the  Cauchy  problem 


i  i 

(•-.  y)  •;;<■•.'/)  Ke|  \\r.l..i\]  Jr  i|  (2.6) 

l**  i  /  :  i  ; 

/  »  i  «* 

satisfying  the  initial  conditions  (see  2.5)  v*(v,0)  -u’(.v)  and  d)  :f‘(v). 


Tn  region  Cl  the  sequence  converges  evenlv  to  an  exact 

i* 


’e solution  of  the  Cauchy  problem  c(x,  y)  for  Laplace 


ii\  -  or.  • 


Now  let  us  construct  an  algorithm  for  tr.c  .-.^merical  solution 
of  fine  Cauchy  problem  which  car.  re  general! z-  c.  a  ea.u  -  or- 


o  a 


complex  equations,  including  nonlinear  ones.  For  this  we  approximate 
the  desired  solution  by  the  Lagrangian  polynomial  for  x: 


u 


/« • 


(2.7) 


where  $jn(y)  —  linear  function  ^(y)  (**.y)  of  values  of  approximate 

solution  for  2n  +  1  rays  of  x  ■  .  Placing  (2.7)  in  (2.1)  and 

requiring  that  the  expression  obtained  is  identically  satisfied 
for  all  rays,  we  obtain  the  approximating  system  of  ordinary 
differential  equations 


GRAPHIC  NOT 
REPRODUCIBLE 


"i; 

*  M 


l) 


(2.8) 


and  initial  conditions 


<»(0) «<»(**).  ^(0)»t(jr,).  (2  9) 

It  is  obvious  that  any  solution  of  initial  equation  (2.1),  which 
at  a  fixed  y  la  a  polynomial  of  tha  power  2n  based  on  x,  for  rays 
x  *  xk  satisfies  the  approximating  system  (2.8).  Since  the 
solution  t*(x,  y)  of  the  Laplace  equation  (2.6)  constitutes  a 
polynomial  for  x  and  y  (for  x  to  the  2n  power),  then 


fU to)  -  *’•  (**•  *)  ** (S  ” 1 5  r+r  1  r‘  +  to)  (2.io) 

/-•  /-* 


constitutes  a  solution  of  the  Cauchy  problem  for  approximating 
system  (2.8)  with  initial  data  (2.9).  This  can  be  easily  checked 
also  by  direct  substitution. 
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Thus,  the  solution  of  approximating  system  (2.8)  with  initial 

data  (2.9)  at  n  -*■  ®  converges  to  an  exact  solution  of  the  Cauchy 
problem  for  the  initial  partial  differential  equation.  Dimensions 
of  convergence  range  are  determined  by  the  distribution  of 
singular  points  for  initial  data  which  are  continued  into  the 
complex  region. 

In  conventional  difference  systems  the  number  of  points  by 
which  derivatives  are  approximated  is  fixed.  Thus,  with  an 
increase  in  the  number  of  steps  the  accuracy  of  the  approximating 
system  increases,  not  due  to  an  increase  of  established  information 
about  the  analytic  properties  of  the  unknown  functions,  but  due  to 
an  improvement  of  the  quality  (accuracy)  of  this  information, 
inasmuch  as  it  is  gathered  from  all  the  smaller  environs  of  the 
examined  point.  In  contrast  to  this  the  difference  system  (2.8), 
a  generalization  of  which  is  used  in  this  work,  may  be  called  an 
"arrangement  with  increasing  information,"  or  an  "analytic" 
system,  since  with  an  increase  in  the  number  of  points  on  a  layer 
there  is  an  increase  of  analytic  information  about  the  function 
which  is  established  in  the  system.  This  most  fully  considers 
the  analytic  nature  of  the  solution.  Let  us  note  that  it  is  not 
obligatory  to  perform  an  approximation  of  derivatives  with  respect 
to  all  points  located  on  the  layer.  It  is  important  that  the 
number  of  these  points  Increases  according  to  a  definite  law 
together  with  a  decrease  of  the  step  in  the  direction  of  the 
Initial  data. 

In  general  the  approximating  system  (2.8)  can  be  integrated 
only  numerically.  During  numerical  integration  the  error  of 
approximate  solution  is  composed  of  three  parts: 

ftj  •(  and  .A*  <r '  -  errors  of  approximation  in  direction  x  and  y 
respectively,  6fe<J>  —  rounding  error.  Let  us  consider  the  properties 
of  growth  of  these  errors.  For  reduction  of  computations  we  will 
consider  simplified  initial  data 

«,  (  r,  0)  --  I(x),  0)  0 

and  will  dispose  N  *  2n  +  1  nodes  on  initial  line  y  *  0  in  zeros 
of  the  Chebyshev  polynomial  ™2n+^(x).  ^he  rernainin£  member  of  the 
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Lagrangian  interpolation  polynomial  in  the  complex  region  c  =  x  +  in 
for  function  a(c)  can  be  represented  in  the  form  of  the  complex 
integral 


6R.V.,„  not 

IWaiCBlT 

where  tL  - 


2:u  J 


P  P» 


T-- 

TV-:  H  »  ). 


(2.11) 


2n+1  -  Chebyshev  polynomial,  p  «  p$  <  p«  -  ellipse,  completely 
disposed  (Pig.  *0  in  region  of  analyticity  of  interpolated  function 
9(0.  Variable  c  is  changed  in  the  region  p<pi<p*-  Parameters 


and 


p^  can  be  taken  as  close  as  desired  to 


Conducting  an 


estimate  of  the  integral  with  the  use  of  the  relationship 


K*‘»  (IV.(*')I)  p. 


according  to  (2.6)  and  (2.4)  we  find  that  for  as  little  a  value  of 
e  >  0  as  desired  there  is  such  an  n*,  that  with  n  >  n* 


I  tip  I  <|/U*)I 


(  ;  *r-« 


.lC.  J 


Here  L2  designates  the  length  of  curve  p  ■  p2,  M2  is  the  naximun  of 
modulus  of  function  o(0  on  curve  p  ■  p2,  and  A  is  the  minimum 
distance  between  curves  p  *  o2  and  p  ■  p^. 


If  It  is  assumed  that  e<  - - (p* — pi), then  according  to  (2.12; 

.S*  -0  at  N  *  2n  ♦  l  ■*•  ».  Thus,  in  the  class  of  analytic  functions 

the  error  connected  with  approximation  in  direction  x  decreases 
exponentially  with  an  increase  in  the  number  of  nodes.  This  Is 
connected  with  the  rapid  decrease  of  error  of  approximate  j  on  of 
initial  data: 


- 1  . 

=  2.1 0, 


2 .  ..  i ) 


Here  r>2  is  the  minimum  distance  of  segment  !  < x  <  1  to 
curve  p  =  p~  and  p^  *>l. 


Now  let  us  consider  the  increase  of  arbitrary  errors  in 
initial  data.  For  this  we  will  take  the  distribution  function  of 
errors  in  mesh  nodes. 


e  whan  X  =*  A, 
0  whao  X  =a 


=*  0, 

I  h  n. 


Initial  data 
following  solution 


«f(.v;.0)-5o(jrj).  f„(xh  0)-0  is  satisfied  by  the 
of  the  approximating  system: 


fl««p  e  Re  (  —  ~~  ^ •  -  *»>  •  •  (*  —  *»-i) (*  (£ -  *>0  \ 

IK  —  *•)  ( r,  -  -  *,) . . .  (jr,  —  *„-»)  (jr„  -  xA  f ,) . . . .  (r„  —  x^)! ' 


Prom  here  it  is  easy  to  obtain  the  asymptotic  expression  of  behavior 
of  error  at  large  n  and  y  >  0: 


2  M 


>ttil 

x  7* +** [x  ws  (2n  t*  !> 0 + y  5,0  (2n  r  I )  01, 


2  -1-  y7» -1  -  par1®. 


(2-14) 


which  corresponds  to  Chebyshev  nodes. 


Thus,  with  an  increase  of  the  number  of  nodes  N  *  2n  +  1 
the  arbitrary  error  in  initial  data  grows  exponentially  with  fixed 
x,  y.  Just  as  when  using  the  usual  "nonanalytic"  difference  systems 
C14]. 


If  during  the  integration  of  the  approximating  system  of 
ordinary  differential  equations  (2.8)  stable  difference  systems 
are  used,  then  increase  of  errors  ;  and  5  is  determined  by 
the  system  of  approximation  in  the  direction  x.  According  to 
(2. 11*)  they  increase  by  exponential  law  with  an  increase  in 
the  number  of  nodes  N  «  2n  +  1.  Error  of  approximation  is 

found  using  a  computer  and  at  fixed  N  can  be  made  sufficiently 
small  so  that  its  Influence  on  the  solution,  taking  growth  Into 
consideration,  does  not  exceed  e. 

Conversely,  the  rounding  off  error  in  each  point  cannot  be 
decreased,  3ince  it  is  determined  by  the  maximum  numbers  of  digit 
positions  utilized  during  calculation,  i.e.,  the  properties  cf  the 
computer.  During  approximation  of  initial  data  the  number  of  nodes 
N  ■  2n  +  1  Is  determined  with  the  required  accuracy  of  approximation 
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of  the  solution  in  this  direction  and  depends  on  the  nature 
of  the  problem.  Prom  this  follows  the  important  conclusion  that 
in  the  elliptic  region  the  Cauchy  problem  can  be  solved  numerically 
only  in  that  case,  when  the  rounding  off  error  is  sufficiently 
small,  so  that  taking  into  account  growth  at  the  selected  N 
they  do  not  exceed  errors  of  approximation. 2  These  considerations 
determine  the  selection  of  parameters  of  the  difference  grid  and 
the  practicable  accuracy  during  solution  of  a  specific  problem. 


S  3 .  Construction  of  an  Algorithm  for  the  Numerical 
Solution  of  a  Boundary  Value  Problem  In  a  MlxecT 
Region,  when  in  One  Part  of  the  Region 
the  Equation  Has  an  fclllptlc  Nature, 
and  the  Other  c  Hyperbolic  Nature 

As  a  model  we  will  consider  the  equation 


ifltp 

<h» 


:  {i 


(3.1) 


Basic  questions  of  the  existence  and  uniqueness  have  been  studied 
for  equations  of  such  a  type  and  more  general  ones  in  the  papers 
[15,  16,  17].  Tf  we  disregard  the  vorticity  of  flow  and  in  the  first 
equation  (1.1)  introduce  the  velocity  potential,  then  (3.1)  can 
serve  as  a  model  of  a  problem  dealing  with  supersonic  flow  past  a 
blunt  body  with  a  detached  shock  wave  in  a  physical  plane.  In 
variables  £,  0,  where 

'c  (O)  rf{  j)  ' 

the  region  of  flow  between  surface  of  body  and  shock  wave  in  Pig. 

1  will  be  converted  into  a  band  (Fig.  5a). 


Let  us  consider  the  band  0  <  x  <  j  (Fig.  5b).  Lines  x  *  0  and 
x  ■  i  correspond  to  shock  and  body  surface.  Regions  ~l<y<l  and 
correspond  to  subsonic  and  supersonic  units  of  flow. 
Coefficients  of  equation  (3.1)  Just  as  the  coefficients  of  the  first 
I'ou.il  Ion  (1.1'  are  oven  functions  of  0.  Equation  of  characteristics 
\  \ .  '  haa  t  ?ic*  same  form 


'J  *  cIj  (jc  — c) 


(3.2) 
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< 

«*.  . 

> 


co  rj» 


Fig.  5a. 


m  it  * 


Fig.  5b. 


Here  c  «  const.  With  x  <  c  we  have  the  characteristic  of  the  family, 
with  x  >  c  —  of  the  secondary.  All  characteristics  touch  the  parabolic 
line  y  *  1,  therefore  the  triangular  region  BCD  (Fig.  5a)  degenerates 
into  a  line  (Fig.  5b).  In  accordance  with  the  formulation  of  the 
gas-dynamic  problem  we  will  consider  as  having  "physical  meaning" 
solutions  with  continuous  and  limited  second  derivatives.  Making 
use  of  the  separation  of  variables  4>  «  F(x)f(y),  we  obtain  two 
families  of  solutions 


const 


•  V ')  K  •'  •'  J  ■  i(  v)  0, 

..  >  •.  r.'Mv) 

f  ■■■<.■)  1  •  o, 


(3.3) 


corresponding  to  plus  and  minus  signs  x£.  From  (3-1)  it  is  clear 
that  a  necessary  condition  of  boundedness  of  second  derivatives  is 
the  requirement  with  y  ■  1.  Considering  the  condition  of 

symmetry  for  axis  at  y  ■  0,  we  obtain  for  solutions  of 

fk(y),  having  "physical  meaning,"  in  the  region  0  £  y  <_  1  cne  boundary 
value  problem 


/'(*)  -0.  MO)  -0. 


(3.^0 


determining  the  fundamental  values  tX^.  For  solution  of  this  boundary 
value  problem  we  use  the  Cauchy  problem  with  initial  data  for  axis 


M0)  '  .  M0)  -o. 


(3.5) 


selecting  by  iteration  e.  such  that  it  satisfies  the  first  condition 
(3.*0. 
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C . 


inot  reproduce 


We  will  start  from  a  case  when  in  (3-3)  the  minus  sign  In  front 

p 

of  A^  is  removed.  From  the  second  equation  we  have 


/’  M;/)  .  >  . 


(3-6) 


and  consequently  |  T  |  1  le^l  on  segment  0  <_  y  <_  1  for  solutions 
satisfying  initial  data  (3.5).  Hence  the  immediate  result  is  that 
the  boundary  value  problem  (3.^)  has  only  a  zero  solution  r,  (v)  z  0 

K 

at  0  <_  y  <_  1.  From  the  general  theory  of  differential  equations 
of  the  type  (3.6)  it  follows  further  that  this  solution,  being 
continued  in  region  y  >  1,  is  also  identically  equal  to  zero. 

Let  us  now  remover  the  plus  sign  in  (3.3)  before  X ,  .  ~n  zr.ij 

/i 

case  the  fundamental  values  of  boundary  value  oroblem  (3.4)  are 

2 

easily  determined  A^  =  k(k  +  1);  as  a  result  we  obtain  a  far.il, v  of 
particular  solutions  of  equation  (3.1)  with  a  limited  second  derivative 


i . ih)  ;U  y  ^*  *6'  *^*‘:hA*  .»'i. 

where  Pf(y)  -  derivative  of  Lagranian  polynomial. 

K 


;  3 . 7 ) 


Let  us  consider  the  question  of  the  uniqueness  o"  t hr-  :ol 
of  boundary  value  problem  for  equation  (5-1).  determined  b\ 
conditions 


■,  ,0 .  v)  •}  (y»  ’j'j  -  'h(y).  v,(fj)  j- 

7his  procler.  constitutes  a  model  of  supersonic  flow  pas 
body  with  a  detached  shock  wave.  Let  us  assure  tr.at  fu; 
6,(y)  and  6^(y)  are  assigned  only  up  to  a  Daraboiic  iir.< 

X  c. 


on  the  segment 


We  will  assume  that  here  tner-e 


solutions  of  equation  (3-1)  with  a  continuous  limited  s 
<b ,  (x,  y)  and  <j> 0 ( x ,  y),  satisfying  boundary  conditions 
their  difference  60  is  solution  (3.1)  with  a  limited  ..c 
derivative,  satisfying  uni  for:,  conditions 


ecorc  ivative 


l 


y 


'NOT  REPRODUCE  ~ 

We  will  look  for  6$  by  solving  the  Cauchy  problem  from.  the  axis 
of  symmetry  y  =  0 

Cf(V.O)  £  (.V),  •  ‘‘  M  ^  '»• 

Here  function  e(x)  is  determined  from  condition  d'bq'dx*  0  at  y  =  1, 
necessary  for  boundedness  of  the  second  derivative  at  y  =  1 . 
Considering  •/. I  *  f.'.V.i)7,  from  the  first  equation  (3.3)  we  obtain 

<\(V)  -f:4w(iUv)  tv). 

Expanding  function  e(x)  in  accordance  with  boundary  conditions 
(3-9)  into  a  Fourier  series  based  on  sines 

e(v) ;  ;  E*  sui  (%.•«). 

we  present  the  desired  solution  in  the  form 

^  }]  /<(y)sia(2*;u). 

*  j 

Functions  f^(y)  should  satisfy  the  second  equation  (3.3) 
and  boundary  conditions  (3.*0.  According  to  the  preceding  it  follows 
immediately  from  here  that  all  f^'y)  =  0  and  in  the  elliptic  region 
0  <_  y  <_  1  =  <t>2  •  If  boundary  conditions  (functions  o,  and  o 

are  continued  up  to  points  G  and  F  (Fig.  5),  then  we  obtain  6.^  =  02) 
in  the  whole  region  AGTFE.  An  analogous  result  can  be  obtained  by 
assignment  on  lines  x  =  0  and  x  =  ^  of  boundary  conditions  of  a  more 
general  form. 

Thus,  ir.  the  class  of  functions  with  continuous  and  limited 
second  derivatives  the  assignment  of  boundary  conditions  (Fig.  5a,  5b) 
on  segments  AB,  ED  and  conditions  of  symmetry  on  segment  AE  uniquely 
determine  the  solution  of  equation  (3*1)  in  elliptic  region  AbDE. 

If  boundary  conditions  are  continued  in  the  nyperbolic  region 
up  to  points  G  and  F,  then  a  unique  solution  is  determined  In  the 
entire  region  AGTFE,  and  consequently  also  in  its  section  AGFE. 

Thus,  the  formulation  of  a  boundary  value  problem  ror  model 
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equation  (3.1)  possesses  all  the  basic  features  of  the  gas-dynamic 
problem  being  examined  (S  1). 

Let  us  consider  the  application  of  the  Cauchy  problem  for 
construction  of  a  boundary  value  problem  (3.8).  Two  primarily 
different  systems  are  possible. 

1.  If  in  addition  to  conditions  (3.8)  on  the  axis  of  symmetry 
y  ■  0  the  distribution  of  unknown  function  $(x,  0)  »  b(x),  is 
assigned,  then  in  band  0<a-<-*-,  y^O  a  mixed  boundary  value  problem  with 
initial  data  at  y  *  0  is  obtained.  Selecting  with  the  help  of 
iterative  process  the  function  b(x)  ■  b^(x)  in  such  a  manner  that 
the  condition  of  boundedness  of  the  second  derivative  is  fulfilled, 
we  obtain  the  desired  solution  of  boundary  value  problem  (3.3). 

This  method  is  close  to  method  of  work  [6], 

We  will  examine  a  deviation  from  an  exact  solution 
•*.7  -,(*•.  y) r*  (.«,//),  caused  by  a  distinction  in  initial  data  of 
Hv)  .  ,  in  the  process  of  iteration.  Function  6<j>(x,  y) 

satisfies  equation  (3.1)  and  the  uniform  boundary  conditions 
(3.9),  so  that  in  accordance  with  the  preceding  it  has  a  limited 
second  derivative  or.]y  at  6$  =  0,  i.e.,  when  the  exact  solution 
is  found.  In  other  words,  since  conditions  (3.8)  uniquely  determine 
the  solution,  then  any  solution  arbitrarily  differing  little  frcr 
the  desired  one  based  on  initial  data  has  an  infinite  second  derivative 
at  y  ■  1.  In  this  method  the  Iterative  process  of  approximation 
to  a  solution  with  a  continuous  and  limited  second  derivative 
is  constructed  in  the  class  of  functions  with  an  infinite  derivative. 
Such  incorrectness  in  the  vicinity  y  *  1  during  the  numerical 
solution  of  a  boundary  value  problem  strongly  aggravates  the  di  "ficul- 
tlns  connected  with  the  usual  incorrectness  of  the  Cauchy  problem 
in  an  elliptic  region.  It  can  be  surmounted  in  a  class  of  sufficiently 
smooth  functions.  As  can  be  seen  from  the  second  equation  (3.3), 
difficulties  increase  progressively  with  an  increase  in  the  narmonic 
number  k,  i.e.,  with  an  increase  of  accuracy  of  the  approximating 
system.  Results  of  this  analysis  are  well  confirmed  during  the 
numerical  solution  of  a  nonlinear  problem  about  supers  or.  i  1  flowing 
around  with  a  detached  shock  wave  [6]. 
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2.  The  second  system  amounts  to  the  problem  being-  examines 
with  initial  data  on  line  x  =  0 


•f(0.  y)  -  'i  j  S' •  y)  -  Mv) 


and  even  function  r.  (>J)  -t?  (y)  will  be  selected  by  iteration  in 
such  a  way  that  the  second  boundary  condition  (3.8)  is  satisfied 
on  line  x  *  1/2.  Since 

y\~\-P\{y)(\-y'). 


>/) 


■  r 


)  (t  -  /) 


BIE  NO' 
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etc.,  then  during  approximation  of  function  i^(y)  by  even  polynomials 
the  approximate  solutions  in  the  process  of  iterations  will  be 
linear  combinations  of  particular  solutions  (3-7).  During  each 
iteration  the  condition  of  symmetry  (third  condition  (3.8))  and 
condition  of  boundeaness  of  second  derivatives  are  satisfied 
automatically,  since  they  are  satisfied  by  the  particular  solutions 

(3.7) .  The  possiblitity  of  constructing  an  iterative  process 

in  a  class  of  functions  with  a  continuous  limited  second  derivative 
in  this  method  is  conditioned  by  the  fact  that  approximate  solutions 
are  not  subordinated  to  any  condition  at  x  =  1/2  (of  three  conditions 

(3.8)  only  two  are  subordinated),  and  consequently,  according 

to  conditions  of  uniqueness,  can  belong  to  this  class  of  functions. 
From  • nis  analysis  it  is  clear  that  during  resolution  of  a  problem 
on  suj  orsonic  flow  past  a  body  with  a  detached  shock  wave  one 
should  use  system  2,  since  it  corresponds  to  the  nature  of  the 
problem. 


Let  us  now  use  the  proposed  difference  system  for  solving 
a  Cauchy  problem  and  a  boundary  value  problem  based  on  system 
?  for  a  model  equation  of  mixed  type  (3.1).  For  this  in  accordance 
with  our  system  we  will  trace  the  straight  lines  y  *  const  =  y^ 
based  on  m  in  the  upper  and  lower  half-planes.  We  will  call  these 
straight  lines  rays.  All  told,  including  the  axis,  we  obtain 
2m  +  1  rays.  Assuming  the  unknown  solution  even  and  aoproxl matinm 
it  with  Lagrangian  polynomials  based  on  y,  we  obtain  for  the  function 

.  ( *.  y)  ~  t*  (*.  y)  :  >]  V;°  i  (  0  y:/ 
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and  for  derivatives  with  respect  to  y 
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(3.11) 


Here  -  linear  function  tmW  !  T**  (t.  y,)  values  of  the 

unknown  function  for  m  ♦  1  rays  in  the  upper  half-plane.  Placing 
(3.11)  in  (3.1)  and  requiring  that  the  resulting  expression  is 
identically  satisfied  on  all  rays,  we  obtain  an  approximating 
system  of  ordinary  differential  equations 


+  <1 -ir)  y  2/(2/-  \) •?%(*) yJl  -  0 

dr 


(3.12) 


(ft  ...  1  a) 

for  determination  of  values  of  approximate  solution  (x,  y) 

on  the  rays. 


One  of  the  evident  properties  of  an  approximating  system 
is  that  if  there  exists  a  solution  of  the  initial  equation  (3.1) 
♦m(x»  y)  which  at  a  fixed  x  is  a  polynomial  of  the  power  2m  with 
respect  to  y,  then  on  the  rays  it  satisfies  the  system  (3.12), 
i.e.,  functions 

‘  '  1  (')  "  •  A  '.i)  v'i  It  ••• 

are  the  solution  of  a  system  (3.12)  with  any  disposition  of  rays. 
This  follows  from  the  fact  that  for  this  solution  •;*,  and 

correlations  (3. 11)  are  exact.  Let  us  consider  a  Cauchy  problem: 
in  region  AGPEA  0  < y^ch(l, -T^)  to  find  the  solution  of 

equation  (3.1),  satisfying  the  initial  data  (see  Fig.  5b): 


(  (0.  y)  f\Jy)  -  a ,  •  j.  . . .  n^y'-\  v,(0.  y)  ). 

Obviously, 

«i/>!(y)(!  -j  )  I*  ...  I- (i/Ml  -y'). 


(3.13) 


where  aQ,  a^,  ...,  are  constants.  The  solution  of  the  Caucny 

problem  for  equation  (3.1)  with  initial  data  (3.13)  will  be  (see  3.7) 
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(3.1*0 


t(x.y)  >J  K> /»1-I (;')<!  -/)clf[/2 i(2j  l)\v)}, 

/I 

and  the  solution  for  the  approximating  system  (3.12)  will  be 

Tta(x)'--*y(jr.fJ  (*-0.1.  ... .  m)  (3.15) 

with  any  situation  of  rays. 

Here  the  solution  of  approximating  system  (3.12)  gives  on  the 
rays  an  exact  solution  of  the  initial  equation  (3.1).  And  if  the 
Lagrangian  interpolation  polynomial  (3.10)  is  used,  then  from  the 
solution  of  the  approximating  system  (3*15)  we  obtain  an  exact 
solution  (3.1*0  in  the  entire  region. 

In  the  case  examined  the  region  of  Influence  of  differential 
equation  (3.1)  is  ARPEA,  while  the  region  of  influence  of  approxi¬ 
mating  system  (3.12)  formally  is  AGFEA ,  since  only  the  values  of 
unknown  function  in  this  region  are  bound  in  the  difference  approxi¬ 
mation. 

It  is  known  that  in  the  hyperbolic  region  in  general  the  region 
of  influence  of  the  approximating  system  must  coincide  or  envelop 
the  region  of  influence  of  the  initial  system  of  differential 
equations.  Otherewise  there  is  no  convergence  of  approximate  solutions 
to  exact,  and  the  difference  system  itself  turns  out  to  be  unstable 
during  calculations.  The  meaning  of  this  affirmation  is  that  if 
the  region  of  influence  of  the  system  of  differential  equations 
exceeds  the  bounds  of  the  region  of  influence  of  the  approximating 
system,  then  by  changing  Initial  data  on  part  of  the  initial  line 
(in  our  case  on  section  GR),  we  do  not  affect  the  solution  of  the 
approximating  system  and  at  the  same  time  we  change  the  solution 
of  the  system  of  differential  equations  In  the  region  of  influence. 

Although  it  is  accurate  in  general,  this  affirmation  becomes 
incorrect  if  it  were  to  remain  in  the  class  of  analytic  solutions. 

This  is  already  clear  from  the  elementary  consideration  that  in 
an  analytic  case  it  is  impossible  to  change  initial  data  on  any 
section  without  having  changed  them  in  the  appropriate  manner  on 
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the  entire  initial  curve,  since  in  this  case  assignment  of  initial 
data  on  any  section  specifies  them  in  the  entire  region  of  analyticity. 
This  position  can  be  given  a  more  thorough  foundation.  Thus,  in  the 
paper  by  Dahlquist  [18]  for  the  wave  equation  during  the  usual 
three-point  system  of  approximation  of  a  derivative  in  direction  of  the 
initial  straight  line  it  is  proven  that  in  a  class  of  analytic 
functions  the  solution  of  an  approximating  system  converges  to  an 
exact  one  also  when  the  region  of  influence  of  a  differential  equation 
exceeds  the  bounds  of  region  of  Influence  of  the  approximating  system 
(i.e.,  at  r/h  >1,  where  h  is  a  step  along  the  initial  data,  and  t 
in  a  perpendicular  direction). 


When  in  a  class  of  analytic  functions  they  are  looking  for  a 
solution  with  the  help  of  an  "analytic"  difference  system,  then, 
and  in  the  hyperbolic  region,  the  regions  of  influence  of  differential 
equations  and  the  difference  system  are  determined  not  by  the  position 
of  extreme  characteristic  and  points  which  are  bound  in  the  system, 
but  by  the  region  of  analyticity  of  the  solution  and  the  convergence 
region  of  the  difference  system.  Let  us  illustrate  this  in  an 
example  of  model  equation  (3.1). 


We  will  construct  the  solution  of  the  Cauchy  problem  with  the 
same  initial  data  (3.13)  in  region  of  influence  of  eauationn  (3.1) 
ARPEA,  bounded  above  by  characteristic  RP,  the  equation  of  which  is 
y  ■  ch(x  -  £q).  In  order  to  use  the  proposed  system  we  will  pass 
from  variables  x,  y  to  variables  x,  n,  where 


ch«-{j  * 


(3.16) 


Here  region  ARFEA  will  turn  into  rectangular 
differential  equation  (3-1)  will  take  the  form 

__  2n  _ **?_.  j _ \  :_j£ _ # .  _ 

Ok*  ch(x  — Oxdi\  <3iJ* 

__  « .  .*?. ,,  o. 

*  cL*  (*-?,)  *1, 


1.  and 


(3.17) 


Just  as  previously,  by  carrying  out,  in  the  upper  half-plane  m  +  1 
rays  n  *  const  ■  (0  £  hj  ±  1)  and  using  approximation  of  the 

unknown  solution  with  even  Lagrangian  polynomials  (3.10),  we 
obtain  approximate  expressions  for  derivatives  with  respect  to  r: 
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in  the  form  (3.11).  Placing  them  in  equation  (3.17)  and  requiring 
that  the  resulting  expression  is  fulfilled  identically  on  the 
rays,  we  obtain  an  approximating  system  for  ordinary  differential 
equations 


• n 

'  *  ch(/-f .)  '  l  /  ‘  j1* 
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(3.18) 


for  determination  of  values  of  approximate  solution  ,  ^f'n  (x, 
y)  on  the  rays.  Since  along  the  extreme  ray  n  =  1 ,  which  is 
a  characteristic  of  differential  equation  0!0x-d!dx.  then  the 
characteristic  equation  for  it  is  written  in  the  form 


v  h:  ( r  !  J  d  <?<f  _  I  •  sh*  (r  -  <'1'.  () 

>h5  “  *.:»  (jr  ■  |yJ  dt  J<\  Ch*(x--|»)  Ji) 


(3.19) 


It  is  easy  to  see  that  the  equation  of  system  (3-16)  for  extreme 
ray  k  =  m,  n  *  1  constitutes  a  difference  recording  of  cnaracterlst i c 
equation  (3.19).  It  is  also  obvious  that  fulfillment  of  (3.19) 
on  line  r  =  1  directly  follows  from  the  differential  equation 
(3.17).  Placing  y  * ijch(x— io),  in  (3-7)  we  obtain  a  system  of  particular 
solutions  of  equation  (3-17).  The  exact  solution  of  cur  problem 
will  obviously  be  (see  3*1*0 

m 

f  (Jr.  y)  «  a,  +  Y4  pU  -i  1*1  ^  (*  —  *o)l  X 

M 

X  1 1  -  .»*  ch  [x  -  J„)J  ch  [  v*i  <2/“T)  *  1} .  ( 3  • 2  J ) 

Since  with  a  fixed  x  this  solution  will  be  a  polynomial  with 
respect  to  n  to  the  Pm  power,  then  the  solution  of  the  approximating 
system  at  any  situation  of  rays  can  be  presented  in  the  iorm 


n{<)  'f(f.'h) 


('<  0,  I ....  ,  in). 


(3.21) 


If  one  were  to  approximate  the  solution  with  the  help  of 
Lagrangian  interpolation  polynomials,  tt°n  from  the  solution  of 
approximating  system  (3.21)  we  obtain  the  c.  ict  solution  of  the 
problem. 
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The  solutions  of  approximating  systems  in  the  form  (3.12)  and 

(3.18)  obviously  agree,  since  they  give  values  of  one  and  the  same 
exact  solution  on  different  lines  (y  ■  const  and  n  ■  const). 

If  one  were  to  assign  a  perturbation  on  initial  line  x  -  0  or 
any  layer  x  -  const,  then  the  growth  of  errors  during  transition 
from  layer  to  layer  and  increase  in  number  of  rays  (growth  m)  will 
be  the  same  for  both  solutions.  In  essence  this  fact  is  obvious 
in  the  class  of  analytic  solutions  (and  consequently  of  initial  data 
also).  Remaining  in  the  region  of  analytic  solutions,  the  example 
cited,  it  is  impossible  to  change  the  initial  data  (3.13)  on  section 
GR  without  changing  them  in  a  corresponding  manner  on  section  AG, 
so  that  the  difference  system  (3.12)  completely  (exactly)  considers 
the  change  of  initial  data  on  the  entire  section  AR. 

In  the  example  considered  due  to  the  presence  of  polynomial 
solutions  of  equation  (3.1)  for  initial  data  (3.13)  it  is  easy  to 
write  out  the  exact  solutions  of  approximating  systems  (3.12)  and 

(3.18) .  However,  the  result  obtained  (independence  of  growth  of 
errors  on  the  selection  of  systems  (3.12)  or  (3.18))  is  not  connected 
with  the  polynomial  character  of  solutions.  For  equation  (3.1)  it 

is  easily  generalized  for  solutions  presented  in  the  form  of  converging 
series  of  polynomial  solutions.  By  a  somewhat  more  complex  path  an 
analogous  result  is  obtained  also  for  model  equations  which  do  not 
have  polynomial  solutions.  This  fact  is  also  confirmed  by  carrying 
out  numerical  calculations  of  flow  past  blunt  bodies  by  the  method 
examined. 

It  is  important  to  note  that  due  to  the  analytic  character  of 
the  difference  system  the  approximating  system  (3.12)  in  elliptic 
and  hyperbolic  regions  has  the  same  properties.  At  a  fixed  value 
x  ■  H  with  an  increase  of  accuracy  of  approximation  with  respect 
to  y  (increase  of  number  of  rays),  errors  in  initial  data  grow 
according  to  the  law  expmH.  During  numerical  determination  of 
the  approximate  solution  of  a  Cauchy  problem  errors  of  approximation 
and  rounding  off  appear.  In  the  examined  case  of  analytic  solution 
and  difference  system  the  errors  of  approximation  on  initial  line 
x  =  0  have  an  analytic  character  anc  decrease  with  a  growth  cf  r 
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according  to  exponential  law,  so  that  at  specific  limitations  on 
location  of  singular  points  of  the  unknown  solution  the  errors 
evoked  by  them  at  x  >  0  do  not  grow  with  an  increase  of  N  *  2m  +  1 
(see  2.12).  Errors  in  rounding  off  have  a  nonanalytic  character, 
therefore  with  an  increase  in  number  or  rays  they  grow  according 
to  exponential  law  (see  2.13).  Since  in  the  class  of  arbitrary 
functions  the  Cauchy  problem  for  elliptic  equations  is  incorrect, 
then  the  exponential  growth  of  rounding  errors  (and  other  arbitrary 
errors)  with  an  increase  in  the  number  of  steps  on  the  initial 
line  cannot  be  eliminated  by  a  change  in  the  difference  system  cr 
the  form  of  the  region  of  influence  of  the  system  in  the  hyperbqlic 
section  of  region  (see  note  on  page  41).  The  number  of  steps  on 
the  initial  line  (number  of  rays)  is  determined  with  the  necessary 
accuracy  of  approximation  of  solution  in  this  direction  and  depends 
on  the  nature  of  the  problem.  In  connection  with  this,  during  the 
solving  of  mixed.  Just  as  of  elliptic,  problems  (see  §  2)  for  the 
construction  of  a  numerical  algorithm  it  is  possible  to  use  the 
Cauchy  problem  only  in  that  case  when  the  rounding  error  is 
sufficiently  small,  so  that  taking  into  account  growth  during  the 
selected  number  of  steps  on  the  initial  line  they  do  not  exceed 
the  errors  of  approximation. 

Considering  everything  said  about  properties  of  a  difl'erence 
system  during  a  discussion  of  formulation  of  the  problem.  One  should 
note  that  use  of  a  characteristic,  limiting  the  region  of  influence 
of  the  system  of  differential  equations  for  gas  dynamics  as  a 
boundary  line  when  writing  the  approximating  system  of  ordinary 
differential  equations  (here  satisfaction  of  the  characteristic 
condition  is  a  result  of  differential  equations,  see,  for  example, 
(3-19),  (3.17)),  does  nothing  either  for  the  formulation  of 
conditions  of  uniqueness  or  from  the  point  of  view  of  improvement 
of  convergence  and  stability  of  the  difference  system  utilized 
here.  On  the  other  hand,  the  transformation  of  coordinates  with 
the  use  of  the  characteristic  leads  to  a  noticeable  complication 
of  the  approximating  system  (see  (3*12)  and  (3.18)),  to  a  dependence 
of  grid  on  M  number  and  other  conditions  of  flowing  around.  Tnis 
hampers  analysis  and  the  use  of  calculation  data. 
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!ow  let  us  rive  an  example  of  the  solution  of  a  boundary  value 
problem  (3.8)  with  coundary  conditions 

.l1-!)  j '  -  Wi  (y),  <p  yj  i'0  \- :  ■  ;  b,-/  -.o  (y)t 

.?,(,o,.o  (3-22) 

for  model  equation  (3.1).  Considering  m  =  2  (3  rays  in  upper  half- 
plane  y  *  0,  y^  and  y^),  we  easily  obtain  the  solution  q*  of  the 
approximating  system  (3.12),  satisfying  the  boundary  conditions 


in  the  form  q#  3  q(x,  y^)  ,  where 


y(*,  y)  =*.*.:  \-Nx 
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(3.2*0 


during  any  situation  of  rays.  In  the  examined  case  the  solution 
of  approximating  system  (3.12)  gives  a  precise  solution  of  initial 
equation  (3.1)  on  rays  q#(x),  and  with  the  use  of  Lagrangian 
interpolation  polynomials  —  in  the  entire  region  q(x,  y)  ,  since  the 
number  of  rays  corresponds  to  the  power  of  polynomials  in  initial 
data  (3.22).  If  initial  data  are  set  by  polynomials  of  tne  2m 
power,  and  number  of  rays  N  in  the  upper  and  lower  half-plane 
(including  axis  y  =  0)  is  less  tnan  2m  +  1,  then  the  solution 
q#(^,  y),  obtained  with  the  help  of  the  approximating  system,  wl  .1 
be  approximate.  This  makes  It  possible  to  study  the  properties 
of  error  6q  =  q  -  q*. 


We  will  limit  ourselves  to  the  simplest  case,  nonetheless 
making  It  possible  to  reveal  the  basic  properties  of  do.  Let  - 
take  rr  =  2  and  examine  again  the  boundary  value  problem.  (2.22), 
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GRAPHIC  NOT 


with  an  exact,  solution  ^ ( x ,  y)  which  is  (3.2*0.  Durinr  detecting 
of  the  approximate  solution  <t>#(x,  y )  we  obtain  the  number  of  rays 
N  =  3  <  2m  +  1 ,  Then,  using  (3.2*1)  it  is  easy  to  obtain 


d'P  AM  I-  jcA.V-1-  P\(y)Vj*  -  J)  [AC^'  1-AbV  -^h,  -t- 
!-  ^Pl(y)  (y*  - 1)  (A  +  A,e-  »'”*). 

AM  a4(l  —  A*). 

A.V  --  2  It4(l  —  A2)  —  AAf], 

1 


A^  =  - 


[«*(-••£> (A  -*»)_ 


->•(:-*■)]• 

AC,  ~ -A')-, AC, 


(3-25) 


where  y  =  h  —  position  of  extreme  ray.  Here  the  error  in  the 
approximation  of  boundary  conditions  is  given  by  expressions 

C'f  (0,y)  -  '!*</■(</  -h'),  A  (3.26) 

Thus,  at  a  fixed  number  of  rays  N  =  3  the  parameter 
q  =  max(a^,  b^)  is  a  measure  of  accuracy  of  the  approximating 
system.  Expression  (3-25)  makes  it  possible  to  trace  the  nature  of 
dependence  of  Sp  on  coordinates,  the  position  of  the  extreme  ray, 
and  accuracy  of  the  approximating  system.  Differentiating  the 
expression  for  <S4>  based  on  h,  will  obtain 


( 

V 


-2Vi  1-  U‘J» 


:  v 


\  '  : 


l)X 


(3.27) 


From  (3.27)  it  is  clear  that  the  weaker  the  dependence  of  6$,  and 
consequently  also  of  the  approximate  solution  4>*(x,  y)  on  the 
position  of  the  extreme  ray  h,  then  the  smaller  is  q  -  max(aJ( ,  b^), 
i.e.,  the  higher  the  accuracy  of  the  approximating  system.  At 
q  '-*■  0  the  solution  of  the  approximating  system  aspires  to  be 
exact,  and  dependence  on  h  vanishes.  From  (3.26)  ann  (3.26)  it 
is  clear  that  6<J>  and  d<f>#/dh  decrease  according  to  linear  law 
together  with  q,  so  that 


(3-23) 


•Vf>  --  *f  («.  j)  -  •  *f*  (*.  y).^  0  (  ^  -) . 

Thus,  |d$*/ah|  (dependence  of  $>*  on  situation  of  rays)  car. 
serve  as  a  measure  of  accuracy  of  the  solution  of  the  approximating 
system. 


§  4 .  Method  of  Solving  a  Nonlinear  Problem 

For  the  solution  of  the  nonlinear  boundary  value  problem 
formulated  in  S  1  we  will  use  the  difference  system  which  was 
studied  in  detail  on  model  equation  in  I  2  and  3.  From  variables 
r,  0  it  is  convenient  to  switch  to  variables  £»  6,  introducing  the 
new  variable 


\ 


7(0) 


-MO)  -MO). 


(4.1) 


rT(0)  and  r 


where ' r 
and  compression  shock. 


r^,(0)  are  the  equations  for  body  contour 


In  these  variables  the  rerion  of  flow 


between  shock  wave  and  body  will  be  converted  into  the  band 
0  <  £  <  1  (see  Fig.  5a) . 


| - 

o 

KM! 


Let  us  construct  for  m  rays  in  the  upper  and  lower  half- 
clar.es,  all  told,  including  the  axis  of  symmetry,  2m  +  1  rays 
(see  Fig.  2). 


For  obtaining  explicit  expressions  we  approximate  the  unknown 
functions  with  Lagrangian  polynomials  and  then  differentiate  them. 
Considering  the  symmetry  of  flow  with  respect  to  the  axis,  we 
obtain  for  the  functions 


■IW.  ■ 

r*  H 
%!•»«*  •> 

CO 


cx 

.  -V 


/  *  /  •  ^ 
and  for  derivatives  with  respect  to  0 

«'  ~>]  2/«J(5)0^  •.  y  r  >](?/  j- /•;.  *  >*  2/i'jv  i. 

‘  ’  /  •.  /  *« 

Values  Uj,  v^,  and  r^  are  linear  functions  of  values  u 

and  r  on  m  ♦  1  rays  in  the  upper  half-plane, 
c 


(4.3) 

v. 


Placing  the  expressions  for  derivatives  with  respect  to  0  in 
(1.1)  and  requiring  that  the  resulting  expressions  are  identically 
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satisfied  on  each  r ay,  we  obtain  the  approximating  system  of 
ordinary  differential  equations  with  respect  to  u^,  v^,  pk» 
of  approximate  values  of  these  functions  on  the  rays.  Solving 
this  system  relative  to  derivatives  with  respect  to  C,  we  finally 
obtain: 

u.j  V 

•  (*u  I-  y*  eta  o4)  -  bko-k  +  («;  -  «,)  + 

is.  .  f  _ +  \  _ 

rfi  *  I. l  *  f,  «/*  J 

“7S?sr+'*,,'*+u*)]i/ 

.pU  '1 _ •»_/• «.  '/i_;  5'*  *»  >  . 

<**1. 
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(4.5) 


‘  :  (';*  1*>0(^*  I’;)'. \  ;• 

‘  '■>;*  !  c-D*; 

subscript  "k"  denotes  values  of  parameters  on  k  ray,  prime::  designate 
derivatives  with  respect  to  0,  enumerable  by  the  formulas  (4.3). 


On  zero  ray  (at  k  ■  0)  all  the  terms  of  the  second  and  fourth 
equations  (4.4)  identically  turn  into  zero.  For  increasing  the 
accuracy  of  the  system  it  is  expedient  here  to  introduce  equation 
for  with  6  =  0.  Differentiating  the  second  equation  for 

0  and  considering  that  with  0=0 
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•  0 


V«0. 


we  obtain 


S>  wj—g  , 

*L  ,w»+s?rj“’ 

'*>  v  V  1  j*  1 M  (<■/•  !•  **■>)• 

>  -l  ■  \>i  L  *'P  J:-o 


(4.6) 


To  system  (it.J()-(4.6)  it  is  necessary  to  add  a  correlation 
connecting  the  departure  of  shock  wave  from  the  body  e  with  angle 
of  inclination  of  the  shock  wave  to  axis  of  symmetry  a: 

•  J J*  «  ”  (*7*  I  $«*)  clg  («*  OJ  -  r'n.  (1|#7) 

Besides  differential  equations,  in  accordance  with  cor  ditons 
of  uniqueness  the  solution  should  be  subordinated  to  boundary 
conditon  (1.3)  on  the  surface  of  the  body  at  £  *  0  and  to  correlation 
for  compression  shock  (1.2)  at  £  *  1.0  (rc  *  r^  +  e)  in  mesh  nodes, 
i.e.,  in  intersection  points  of  the  ray  with  surface  of  the  body 
and  the  shock  wave.  The  algorithm  for  the  numerical  solution  of 
the  problem  amounts  to  the  following.  Assignment  (m  +  1)  of 
parameters  r^  (J  ■  0,  1,  ...,  m)  according  to  (4.2)  determines 
approximately  the  equation  of  the  shock  wave,  and  with  the  help  of 
correlations  (1.2)  all  the  gas-dynamic  parameters  behind  the 
compression  shock.  Then  by  solving  the  Cauchy  problem  for  the 
system  (lj.JO-(4.7),  we  determine  the  values  of  parameters  in 
nodes  on  the  surface  of  the  body,  which  generally  speaking  do  not 
satisfy  boundary  condition  (1.3).  Selecting  by  iterations  the 
values  r?  in  such  a  way  that  in  all  the  nodes  on  the  surface  of  the 
body  the  boundary  condition  (1.3)  is  fulfilled  with  the  assigned 
accuracy,  we  obtain  the  unknown  solution  of  the  approximating 
system  in  m  approximation. 


Thus  the  method  of  calculation  amounts  to  the  following, 
construct  a  sequence  of  solutions  satisfying  all  the  boundary 
conditions  in  region  AGPF. :  correlations  for  compression  shock 
AG,  conditions  of  symmetry  on  a  segment  of  axis  AK,  and  the 
boundary  condition  on  contour  of  body  AF. 

An  additional  condition  of  boundedness  of  derivatives 
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'nunconversicn  into  zero  of 


Jacobian 


j  (p„  ty) 


is  fui: 


,*ec 


automati¬ 


cally,  since  each  term  in  the  sequence  of  approximate  solutions  satis¬ 
fies  it.  If  at  the  beginning  there  is  assigned  or  if  in  the  process  of 
iterations  there  develops  a  form  of  shock  in  which  in  the  supersonic 
segment  of  region  AGFE  a  limiting  line  emerges,  then  it  will  auto¬ 
matically  be  related  by  the  program  since  the  process  of  calculations 
is  disrupted.  The  desired  solution  and  the  one  obtained  in  the 
process  of  iterations  belong  to  the  class  of  analytic  functions 
where  the  Cauchy  problem  is  correct. 


S  5 .  Examples 

Without  dwelling  here  in  detail  on  results  of  investigating 
the  attributes  of  the  method,  growth  of  errors,  and  convergence 
and  stability  during  the  solving  of  a  nonlinear  problem  on  super¬ 
sonic  flow  around  with  a  detached  shock  wave,  we  will  give  examples 
illustrating  the  solution  of  different  problems  of  flow  around  with 
a  detached  shock  wave. 


Figure  6-9  Illustrate  the  influence  of  M  number  and  body 
form  on  the  geometric  picture  (relative  position  body  surfaces, 
shock  wave,  and  sonic  line)  of  flow  around  by  an  ideal  gas  with 
Y  =  1.4.  Figure  6  shows  the  Influence  of  a  change  in  M  number  over  a 
wide  range  during  flow  past  on  ellipsoid  with  a  ration  of  semi- 
axis  6  *  2.0.  Figure  7  gives  a  geometric  picture  of  flow  at  M  ■  3 
past  a  family  of  bodies  with  front  section  contours  assigned  by 

the  equation  xn  ♦  yn  ■  1.  At  n  >  2  radius  of  curvature  of  contour 

*  >■ 

In  a  critical  point  RQ  is  equal  to  infinity,  and  with  a  change  of 

n  from  2  to  10  the  ratio  of  the  minimum  angle  in  the  vicinity,  the 

H 

radius  of  curvature,  to  body  diameter  ^ln  is  changed  from  0.5 

to  0.07.  Figure  8,  in  a  example  of  flow  at  M  ■  3  past  Cassinian 

2  2  2  2  2  2 

ovals  with  the  contour  equation  (x  ♦  y  )  ♦  2c  (x  -  y  )  * 

4  4  2  2 

*  a  -c(a  ♦  c  -1),  shows  the  influence  of  concavity  of  contour 
in  region  of  critical  point.  Figure  9  shows  flow  at  M  »  3  past  a, 
the  contour  of  which  har  a  discontinuity  of  curvature  at  point  C 
and  is  formed  by  conjugate  circumferences  of  radii  Rq/D  =  1  in  region 
of  critical  point  R-^/D  ■  0.2  in  the  vicinity  of  the  midsection,  and 
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past,  body  F,  constituting  a  ;•  3-dc  ~re-o  segment  with  ar.  ar.rular  point. 
Figure  10  fives  the  distribution  of  sressure  or.  tr.e  surface  of 
bodies  of  various  for r.  at  flow  past  with  '•'  =  3-  f o  rhat  the  curves 
do  rot  merge,  polar  anrle  6  is  accepted  as  an  argument  along  a  he 
surface  of  the  body  (see  Fig.  6,  7,  3).  For  body  f  £  =  6..  (where 
0  --  angle  with  pole  in  center  of  ci rcur.ference  of  radius  r"  )  . 

"or  body  a  up  to  a  point  with  6  =  G  ,  and  further  9  9.,  + 

(see  Fig.  9).  Points  plot  the  experimental  data  of  Yu.  Ya.  oarr  oysl:  i  v . 


Fig.  3.  Fig.  9. 


1  r.vest igation  of  supersonic  flow  past  a  spnere  'ey  r.cnecuili- 
n I urn  air  was  performed  taking  into  account  kinetics  cf  excitation 
o fluctuations  in  and  0^,  reactions  of  dissociation  ana  exchange, 
including  0o  +  =  2N0),  and  ionization  with  the  help  of  atomic 

collisions  N  +  0  =  N0+  +e,  N  +  N  =  N*  +  e,  0  +  0  =  0*  +  e , 

C  c~ 

satisfactorily  describing  the  process  up  to  M  'v  30-  Influence  of 
physical-chemical  transformations  occurring  behind  the  shock  waves 
on  the  departure  of  shock  wave  is  shown  in  Fig.  11. 


Fig.  11. 


On  Figs. 12  and  13  are  shown  the  geometric  picture  and  distribution 
of  pressure  on  the  surface  during  transverse  flow  past  a  cylinder 
(two-dimensional  problem)  with  different  M  numbers. 


(NOT  REPRODUCIB 


Tn  work  [SI"'  a  r.ethod  is  Generalized  f  tm  c  re-  of  :,raf ial 

flew  past  blunted  bodies  with  detached  sh.o  k  wave .  'Inure  It 
shows  the  influence  of  anirle  cf  attack  on  the  re cret.  ic  pic*. are  ef 
flow  past  an  ellipsoid  with  6  =  3.0?  at  V  =  3.  For  tne  same  car 
r’’lg.  15  shows  distribution  of  pressure  on  the  surface  of  *•  bony  In 
a  plane  of  svmretry.  The  points  are  plots  of  the  experimental  data 
of  vu.  Ya.  ’darpeyskiy. 


Pip.  13.  Pip.  1**. 


The  calculations  presented  in  Figs.  6-10  and  IP  and  13  w'-r-o 
conducted  on  a  d-polnt  systen  (selection  of  5  parameters )  ,  nr.u  in 
Fir.  11  on  a  5-point  system  (selection  of  3  parameters ) .  burin' 
calculations  of  spatial  flow  past  during  iteration.:  nf  tne  *'or 
compression  shock  the  selection  13  parameter.'  was.  u  r^ir-f  -..  t 
out  swell  isr  .acre  on  tne  results  of  invest:  at!  a"  r.':,'  , 


Fin.  15. 


note  that  total  error  during  calculation  of  the  majority  of  variants 
does  not  exceed  1% .  For  example,  in  the  case  of  flow  past  a  sphere 
•  V  an  ideal  ras  at  N!  =  3  the  total  error  already  durinr  calculation 
based  on  a  5-point  system  does  not  exceed  0.2f  in  the  entire  field 
1  etween  wave  and  body.  During  calculation  of  flow  past  a  body  with 
a  :’racture  of  the  generatrix  in  the  sonic  point  (B  in  Fir.  9)  '  n 
the  vicinity  of  singular  point  an  asymptotic  solution  was  used. 

"•'his  was  obtained  by  Guderley  [1^]  for  flat  potential  flows . 

hen  the  renerali  zatlon  in  a  number  of  works  [21,  20]  for  vortex 
axi symmetrical  flows  was  used. 


The  authors  thank  G.  I.  Petrov  and  L.  A.  Chudov  for  useful 
discussions  and  consideration  of  results. 
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Footnotes 


'Here  two  cases  of  disposition  of  nodes  are  examined.  For 
an  example,  of  theorems  of  convergence  in  a  general  case,  see 
[12,  13]. 

2  At  present  L.  0.  Frank  and  L.  A.  Chudov  are  developing  methods 
based  on  the  smoothing  of  the  solution  for  each  aye r,  thus  making 
It  possible  to  lessen  the  growth  of  rounding,  errors  [22]. 
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